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Abstract 

Magnetic bearings are subject to performance limits which are quite different from 
those of conventional bearings. These are due in part to the inherent nonlinearity 
of the device and in part to its electrical nature. Three important nonideal beha- 
viors are presented: peak force capacity, force slew rate limitation, and sensitivity 
to rotor motion at large displacements. The problem of identifying the dynamic 
requirements of a magnetic bearing when used to support a known structure subject to 
known loads is discussed in the context of these limitations. Several simple design 
tools result from this investigation. 


Introduction 

Magnetic bearings are moving from the realm of science fiction into that of practical 
engineering. Not only is their feasibility being widely demonstrated, but many 
advantages over conventional bearings are becoming apparent. Amidst this atmosphere 
of optimism, it is important to recognize and understand the shortcomings of this 
class of devices. Hopefully, if our understanding of the limitations of magnetic 
bearings can at least keep pace with our concept of their advantages, premature 
disappointment can be avoided. 

To the engineer accustomed to conventional bearings, the limiting behavior of 
magnetic bearings will be very unfamiliar. In stark contrast to fluid film bearings, 
these bearings become softer (less stiff) as the shaft excursion approaches the 
bearing clearance. Additionally, and perhaps even more alien, the rate at which the 
bearing force can change is strictly limited. This phenomenon is referred to as 
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force slew rate limitation. Like the force slew rate, the peak force which can be 
generated is also subject to strict limits. These characteristics, which represent 
nonlinearities, are especially insistent because their onset can be quite abrupt. 

Very little attention has been focused on these limitations in the literature. 
Britcher (ref. 1) and Sarma and Yamamura (ref. 2) have discussed operation of bear- 
ings in the magnetically saturated regime, for applications where field strength 
requirements are sufficiently stringent to justify the complexity of compensating for 
the resulting nonlinearity. Hebbale (ref. 3) presents a fairly thorough investiga- 
tion of nonlinear performance in the unsaturated regime, but the study concentrates 
on eddy current effects in unlaminated rotors. Lamination of the rotor has become 
common practice and the eddy currents are readily reduced to where their contribution 
to nonlinear effects is minimal. 

This paper represents an effort to describe some aspects of the limiting behavior of 
magnetic bearings and to suggest methods by which the required dynamic capacity of 
the bearing may be estimated. 


Nomenc latur e 


A area of each air gap perpendicular 
^ to the magnetic flux 

A minimum pole cross sectional area 
^ perpendicular to the magnetic flux 

B magnetic flux density 

B saturation flux density 

sat 

bearing force in the x direction 

F bearing load capacity 

max 

f force slew rate: dF/dt 

g length of the air gap 

G nominal air gap: g(x = 0) 

I, coil bias current: invariant 

b 

i coil perturbation current 

magnetic actuator gain: dF^/di^ 

K nominal bearing stiffness 

nom 

k feedback gain: i /x 

o P 

open loop bearing stiffness 

closed loop bearing stiffness 


L coil electrical inductance 

c 

N number of coil turns per pole 

pair 

R coil electrical resistance 

c 

Rj. current sensing resistance 

t time 


voltage across coil 

V g power supply voltage 

x displacement of the rotor from 

the centered position 
x^ geometric clearance 

P sensitivity of the air gap to 

shaft displacement: dg/dx 
nr nondimensional bearing parameter 

6 bearing characteristic length 

magnetic permeability of air 

a nondimensional displacement: x/6 

<f> magnetic flux 


Bearing Description 

A wide variety of magnetic configurations has been proposed for accomplishing magne- 
tic levitation. In applications where load capacity must be maximized relative to 
the size of the bearing package, the configuration described by figures 1 and 2 is 
the most commonly used. (refs. 3,4,5 f 6) For such a bearing, the force generated is a 
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( 1 ) 


nonlinear function of the current in the coils and the rotor displacement: 

F x = ^ 0 P N2A ( (ii/gi) 2 - (iz^) 2 ] 

where the gap lengths have been assumed to vary linearly with rotor displacement- 

gj = G - /3x, g 2 = G + Px (2) 

The bearing is made to do useful work on a system by varying the currents i t and i 2 
in a manner correlated to the rotor displacement. Various methods have been proposed 
for linearizing the relation between a control signal and the resulting bearing 
force. These include establishing a bias flux and using the control signal to per- 
turb it (refs. 5 and 6). relating i* and i 2 to the square root of the control signal 
(ref. 5), and using feedback of the actual magnetic flux level (ref. 4). In the 
discussion which ensues, the bearing is assumed to be bias flux linearized. The 
resulting performance and nonlinearities are similar to what would be obtained with 
square root linearization. Flux feedback has the potential to greatly reduce flux 
saturation nonlinearities and eliminate softening of the bearing with large displace- 
ments. However, no actual applications of this linearization have been reported; it 
seems more useful to discuss the prevalent designs. 


Two methods have been reported for accomplishing bias flux linearization. The more 
common of the two is to establish a bias current in each of the bearing coils and 
then modulate i t by adding a perturbation current proportional to the control signal. 
The other scheme is functionally identical except that the bias flux is established 
by linking the flux from a set of permanent magnets to that of the bearing coils in 
the air gaps. The latter scheme is somewhat more complicated to construct but can be 
substantially more energy efficient. Either scheme, however, is governed by the same 
operating equation and is subject to the same limiting behavior. 


If bias flux 


linearization is employed, then equation (1) 


F = 
x 


u A N 2 /3 

_2_S 


( I K + 1 V 

v b p' 

(G - px) : 


(I,- 1 ) 

v b p' 
(G + Px) 


can be recast as 
2 ' 

2 


(3) 


As mentioned above, permanent magnets can be used to establish the bias flux rather 
them bias currents. In this case 1^ would be replaced by the field strength of the 

permanent magnets, scaled appropriately. 


Maximum Bearing Force Limits 


The most immediately apparent limiting behavior of magnetic bearings is that of peak 
force limitation. This limitation is primarily due to nonlinearity in the magnetiza- 
tion curve of the electromagnet core material. Equation (1) is based on the assump- 
tion that the core magnetic flux is simply proportional to the magnetomotive force 
(MMF: coil turns x coil current). The actual relationship between these two quanti- 
ties is depicted in the magnetization curve for the core material. Figure 3 shows a 
typical curve for conventional magnet iron. The significant feature is the leveling 
of the curve at high MMF’s; beyond a certain point, the flux density ceases to in- 
crease substantially. This phenomenon is referred to as magnetic saturation and the 
flux density at which it occurs is the saturation flux density. Typically, this 
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number ranges between 1.2 and 1.6 Tesla, while some special materials push this limit 
slightly over 2.0 Tesla. The effect of magnetic saturation is that, once a certain 
coil current is reached, further increases in current will produce relatively little 
increase in bearing force. 


The peak force which cam be developed by a given magnetic bearing is determined in 
the following manner. The maximum flux in the most restricted cross section of the 
pole structure is 


<t> 

P 


= B A 
sat p 


(4) 


The maximum flux in each air gap of the magnetic circuit is 


$ 

g 


p NI A / 2G 
o sat g 


(5) 


Since these two fluxes must 


Combining equations (6) and 


be equal, equations (4) and 


NI = 2GB A / p A 
sat sat p o g 


(3) by setting i = 1^ , I^= 


(5) can be combined to obtain 

( 6 ) 

l AI , and x = 0 yields 
sat J 


F = 
max 


PB 2 A 2 /(i A 
sat p o g 


(7) 


It is readily shown that F is maximized by setting I, = l AI . Throughout the 
J max J & b sat & 

remainder of this paper, it will be assumed that this design rule is followed. 


As an example, consider a bearing designed with a pole face area of 10 cm 2 (1.6 in 2 ), 
a minimum pole cross sectional area of 7 cm 2 (1.1 in 2 ), and a geometry factor of 
0.93. If the saturation flux density is 1.2 Tesla, equation (7) indicates that the 
peak force which can be developed is 524 N (118 lb). 


Force Slew Rate Limits 


The bearing force slew rate is limited because the magnet coils have a high induc- 
tance and the power supply to the driving amplifier is at a fixed voltage. For small 
displacements about x = 0, equation (3) can be linearized as 


where 



K x 
x 


( 8 ) 


K.= p A N 2 I, P / G 2 and K = - p A H 2 I?p 2 / G 3 (9) 

log b x o g b v J 

If x is held constant so that we can investigate how rapidly the bearing force can be 
changed by varying the current then 


dF /dt I ^ = K, di /dt 
x *x=0 i p 


( 10 ) 


The voltage across the electromagnet coil is 



If the output stage of the amplifier which drives the coil Is described by figure 4, 
where it is essentially a variable resistance between the coil and ground (the resis- 
tance being variable between about 0.5 ohms and nearly infinity), then the fastest 
possible positive force slew rate is 

dF/dt < K. [ V s - (I b+ i p )(R c+ 0.5 + R f ) ]/L c (12) 

For this output stage configuration, the largest negative force slew rate is 
determined by the reverse breakdown voltage of the output transistor. This 
limitation is typically larger (in magnitude) than that imposed by equation (12). 
Similar arguments apply to P-channel or bipolar output devices, push-pull 
configurations, and pulse-width-modulated amplifiers; the slew rate limitation Is 
inherent in the fixed power supply voltage V . 

Figure 5 illustrates this effect. The trace shown is of the coil current actually 
delivered to the bearing coils when the slew rate of the sinusiodal control signal 
exceeds the maximum slew rate permitted by the power supply voltage (200 VDC) and the 
coil inductance (0.95 H) . Notice that the wave form is distorted only on the 
ascending portion. Distortion begins in the region where the positive slew rate is 
greatest, but is more evident as time proceeds because the current error is the time 
integral of the slew rate error. Once the slew rate limitation sets in, the current 
simply increases at the rate set by equation (12) until the demanded current and 
actual current once again coincide. In the case illustrated by figure 5, this occurs 
slightly after the peak. The descending portion of the waveform is undistorted 
because the reverse breakdown voltage of the output transistor (900 V) is far in 
excess of the 200 V power supply. 

Slew rate limitation has two effects on the performance of the bearing. First, since 
it causes the bearing force to change more slowly than the control signal demands, it 
introduces phase lag. This is shown in figure 6. where the phase shift through the 
amplifier is plotted as a function of frequency for various power supply voltages. 
This produces a sudden loss of damping at high amplitudes, which can be catastrophic. 
The other effect is due to the asymmetry of the distortion relative to the bias 
level. This effectively reduces the bias current which, in accordance with equation 
(9), reduces the actuator gain. K . The result is a reduction in both stiffness and 
damping of the bearing. 

In order to understand how to design around the force slew rate limitation, equation 
(12) must be examined with an eye to what the significant design parameters are. The 
inductance of the magnet coil is defined by 

L = (N*)/I (13) 

c 


where I is the coil current required to produce a magnetic flux <f>. Equations (4), 
(5). (7). and (13) can be used to convert equation (12) to the form: 


dF /dt < 
x 


2j3B A 
K sat p 

p A N 
o g 


V - (I K + i )(R + 0.5 + 


V ] 


(14) 


In general, the term (I^+ lp)(R c + 0-5 + R^) is on the order of 5 A of V^. The ratio 

A /A should be 1.0 for a thrust bearing (0 = 1) and about 0.99 for radial bearings 
P S 
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(/} ~ 0.93). Thus the controlling design parameters are the power supply voltage and 
the number of coil turns. Requiring a large power supply voltage will degrade the 
thermal efficiency of the power amplifier which drives the bearing. However, space 
restrictions and thermal requirements will generally limit the usable number of coil 
windings. Both parameters must be carefully juggled in order to obtain adequate slew 
rate capacity. 


Displacement Sensitivity 


Displacement sensitivity of magnetic bearings is most evident as a softening of the 
bearing, or reduction of its spring rate, as the shaft approaches the radial clear- 
ance. This is a consequence of the nonlinear nature of equation (3). In order to 
readily describe this behavior, it will be assumed that the bearing has been designed 
to operate as a simple spring, without damping. To accomplish this, the perturbation 
current, i , is controlled by feedback of the shaft position: 


The force equation (3) becomes 


I = -k x 
p o 


(15) 


V A N p 
F = -2-fi 

X 4 


(I, - k x) 2 (I,+ k x) 2 

v b o * _ v b o _ 

L (G - Px) 2 (G + 0x) 2 


(16) 


The effective spring rate is the negative derivative of the bearing force with 
respect to deflection: 


p A N p 

K - - dF /dx = - 2 - s (k G - 01, ) 

x,cl x 2 v o ^ b y 


I, - k x I, + k x 
b o b o 


(G - Px)' 


(G + Px)' 


(17) 


Equation (17) is more tractable when nondimens ional ized . To do this, two bearing 
characteristics are introduced: the maximum force or load capacity, F , and the 


max 


nominal stiffness, K 


nom 


Noting that the nominal stiffness is that where x = 0, 


K = |i A N 2 BI U (k G - )3L ) / G c 
nom o g b v o b J 


(18) 


If the bias current is half of the saturation current then the maximum force will be 

F =M N 2 /3Lk /(k G + j3I u ) 2 (19) 

max og bo o b 


generated when x = -I, /k : 

bo 


Note that equation (19) will predict a lower maximum force than will equation (7), 
where the maximum force was found at x = 0 and without defining a feedback law. 
Equations (18) and (19) define a characteristic length of the bearing 


5 = F / K 

max nom 


( 20 ) 


If the following nondi mens ional quantities are defined: 


7 = J3L / k G , r(nr) = (1 +t) (l-nrVnr , a = x/6 , k = K . / K 

bo k j \ j \ j x.cl nom 


( 21 ) 
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then equation (17) can be nondimensionalized as 


K 



nrr - a 


(r - a ) 3 


+ 


nrr + a 
(r + a) 3 . 


( 22 ) 


The definitions of equation (21) lead to a useful expression for the air gap in terms 
of the characteristic length, 6, and the bearing parameter, *7* 

G = (36 ( l+'r) z ( I- tJ/t (23) 

Equation (22) is particularly useful because it describes the variation of bearing 
stiffness in terms of the nondimensional shaft displacement, a, and the bearing para- 
meter, t, in such a way that a and t are entirely independent. This permits compari- 
son of a wide variety of bearings having the same stiffness and load capacity but 
different clearances. 


Figure 7 illustrates the family of curves defined by equation (22) for 0 < nr i 0.66, 

0 < a i 1.0. Stiffness curves for ~r = .64, .50, .40. and .10 are shown separately in 
figure 8. The geometric requirement that G £ 6 dictates [by equation (23)] that 
0 < *r < 0.755. Further, since magnetic saturation sets in for |fi| > 1, equation (22) 
is not valid for values of 6 outside this range. Clearly, the bearing is most linear 
in the region of t = 0.5 and as tt approaches zero. In interpreting these graphs, it 
is important to recognize that equation (22) implies that, as tt approaches zero, the 
air gap, G, becomes very large. This means that establishing the required magnetic 
flux in the air gap will require very high coil currents with the accompanying high 
I 2 R power losses. Leakage flux losses are also accentuated by large air gaps. 

To illustrate the use of these equations consider the design of a thrust bearing 
(P = 1.0) which is to have a nominal stiffness of 2xl0 6 N/m (l.lxlO 4 lb/in) and a 
load capacity of 500 N (112 lb). The characteristic length, 6, is 0.25 mm £0.010 in). 
If a bearing parameter, 'r, of 0.5 is selected for high linearity and low I R losses, 
then the air gap would be [from equation (23)] 0.56 mm (0.022 in). At a rotor dis- 
placement of 0.20 mm (0.008 in), equation (22) predicts that the stiffness would be 
88% of the nominal stiffness, or 1.8xl0 6 N/m (9700 lb/in). At 0.25 mm, this figure 
would decrease to 64%. 


Estimating Dynamic Requirements 

The preceding discussion illustrates the importance of staying within the performance 
bounds of the magnetic bearing. Exceeding the peak force or slew rate limitations 
can introduce a sudden drop in bearing stiffness and damping. The resulting dynamic 
behavior can be disastrous. If the bearing is designed to minimize its displacement 
sensitivity, then it may be permissible to neglect this additional nonlinearity, 
especially if the design load capacity is somewhat greater than the actual dynamic 
requirement of the bearing. For any of these considerations, a careful assessment 
must be made of the actual forces that the bearing will be asked to deliver. 

In estimating the dynamic requirements of a magnetic bearing, the designer may initi- 
ally incline to set them equal to the worst dynamic characteristics of the antici- 
pated loads. However, this is unrealistic and commonly demands far more performance 
of the bearing than is needed. In the following, two types of loads are discussed: 
step and sinusoidal. For the sake of simplicity, it will be assumed that the system 
being supported by the bearing is a simple mass. This assumption permits some fairly 
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general conclusions to be reached concerning slew rate and peak force requirements 
which will provide useful first pass design information and, more importantly, will 
help develop a sense of where these requirements come from. More exact estimates of 
the dynamic requirements placed on the bearing require detailed knowledge of the 
supported structure and the control algorithm. 

Assume that a rigid body of mass m is supported by a control 

geometric restrictions, the body cannot be permitted to move 
from its unloaded position. If a step force of F is applied 
time trajectory is described by 

m d 2 x/dt 2 = F^- F g : x(0) = x q , dx/dt(0) = 0 (24) 


force, F . Because of 
c 

more than a distance x 
to this body, then its 


If, further, the control force is assumed to be initially zero and then changes at a 

fixed slew rate of f- F^= ft, then equation (22) can be solved as 

x = ft 3 /6m - F t 2 /2m + x , t > 0 (25) 

so" v ' 

The body will collide with the geometric constraint if x = 0 for some t > 0. It can 
readily be shown that equation (25) has no positive real roots in t if 

f > (2FV3mx o ) 14 (26) 

For such a simple control scheme (essentially the integral of bang-bang), the peak 
force required of the controller is twice the applied load. If the force slew rate 
were unrestricted, then the peak force need be only slightly greater than F . 

As an example, consider a body having a mass of 0.5 Kg (1.1 lbm) . If the largest 
allowable excursion from its undisturbed position is 0.5 mm (0.020") and a step load 
of 100 N (22.5 lb) is applied, then a force slew rate of at least 52,000 N/sec 

(11,600 lb/sec) is required to prevent a collision. Figure 9 shows the time response 

of this system. The control algorithm used in this simulation consisted of 
delivering a constant force slew rate of ±52,000 N/sec or zero. 


When the same system is excited by a sinusoidal force at frequency u, then the analy- 
sis is even simpler. The best control force will obviously be sinusoidal, 180° out 
of phase with the load. The time response of the system is 


x = (F c ~ F s )/(ma> 2 ) sin<Jt 


(27) 


The restriction that |x| < x^ implies that 

F - x mo) 2 < F < F + x moj 2 
so c s o 


(28) 


For a sinusoidal control force, the peak force slew rate is simply u times 
forced 

wF - x mw 3 < f < (i)F + x mu 3 
s o so 


the peak 


(29) 
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If the same mass as in the preceding example were acted upon by a force of constant 
amplitude 50 N (11.2 lb) but over a frequency range of 100 rad/sec to 1000 rad/sec 
(955 RPM to 9550 RPM) then the control force actuator would have to be capable of 
delivering a force as large as 47.5 N (10.7 lb). The maximum required slew rate 
would occur at 

a(o>F s - x o t™ )/3u = 0 => w = 258,2 rad/sec 

Thus, the actuator must be able to provide a slew rate of 8,607 N/sec (1,934 lb/sec). 


Conclusion 

This paper has treated three forms of nonideal behavior found in magnetic bearings. 
These include peak force limitation due to magnetic saturation, force slew rate 
limits due to magnet coil inductance and finite power supply voltages, and displace- 
ment sensitivity due to the nonlinear dependence of the bearing force on rotor posi 
tion. The mechanisms which give rise to these limitations have been explored in some 
detail with the purpose of developing several relatively simple equations which can 
be employed in the design process. 

As a complement and, perhaps, to motivate this interest in bearing limitations, a 
very simple discussion was given of the minimum dynamic requirements made on a bear 
ing by a sinusoidally or step excited mass. The intent here was twofold. First, it 
was desired to demonstrate that the dynamic requirements do not necesarily match the 
characteristics of the external loads. Second, this simple analysis provides some 
useful guidelines for estimating the requirements of actual systems. 

When working with a mechanism as readily tailored as a magnetic bearing, it is easy 
to lose sight of its limitations. These devices present the engineer wi th a solution 
to a vast range of difficult bearing problems, but if misapplied, their performance 
can be very disappointing. Hopefully, this investigation will help provide methods 
for ensuring that magnetic bearing designs will perform as intended. 
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